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ABSTRACT 

In this paper we present and discuss the effects of scattered light echoes (LE) on the 
luminosity and spectral appearance of Type la Supernovae (SNe). After introducing 
the basic concepts of LE spectral synthesis, by means of LE models and real observa¬ 
tions we investigate the deviations from pure SN spectra, light and colour curves, the 
signatures that witness the presence of a LE and the possible inferences on the ex¬ 
tinction law. The effects on the photometric parameters and spectral features are also 
discussed. In particular, for the case of circumstellar dust, LEs are found to introduce 
an apparent relation between the post-maximum decline rate and the absolute lumi¬ 
nosity which is most likely going to affect the well known Pskowski-Phillips relation. 

Key words: supernovae: general - supernovae: individual (SN 1991T) - supernovae: 
individual (SN 1998bu) - supernovae: individual (SN 2001el) - ISM: dust, extinction 
nebulae 


1 INTRODUCTION 

Due to the crucial role played by Type la SNe in modern c os- 
mology (see for example the review bv iLeibundgutlEoOlIl , it 
is essential to understand the systematics which might affect 
the recent conclusions about the expansion rate of the Uni¬ 
verse. The availability of larger telescopes and the renewed 
interest for these objects as cosmological probes has opened 
the possibility of studying in great detail their evolution, an 
investigation which in the past was feasible for the closest 
events only. In these last few years, in fact, numerous la have 
been followed with an extensive time coverage as late as one 
year past maxim um light, leading to very interesting results 

(see for example iBenetti et al.ll200f^^ ■ _ 

I n two cases, namely those of SNM991T_JSchnndte^]J 
|l994i) and SN 1998bu iCappellaroe^^ 
ISpvromilio et al .l200.dl . the spectra at advanced phases were 
clearly dominated by a phenomenon which had nothing 
to do with the explosion mechanism but rather with the 
neighboring environment. In both cases this very atypical 
behaviour could be explained with the contamination by 
SN light scattered into the line of sight by dust present 
in the surroundings of the explosion. The presence of 
such a phenomenon, also known as light echo (hereafter 
LE), was later confirmed by high spatial resoluti on HST 
imaging llSparks et al.lll9^ : iGarnavich et al.ll20nill . Other 
cases have been discovered for core-collaps e SNe like 
SN 1987 A ICrotts, Kunkel &: McCarthyl [l98^, SN 1993J 

iSueerman fc Crott j l2002l : iLiu. Bremiari fc SeitzeJ l2003ll . 


SN 1999ev^_jMaun^^^^n^^_|200^)__a^ SN 2003gd 
iSugermanI 1200,’ll : IVair^vkTrfcMlippenkol boOfjl . In the 
most recent years, L Es have been searched a l so in the case 
of h i storical SN e jRo maniello et alJ l200!Tt iKrause et alJ 
1200111:^,681 et al ]|20 oT 

With the purpose of studying this effect and its possible 
applications to the analysis of the SN environments and their 
relation with the SN progenitors, we h ave started a series of 
papers on this topic. The first of them (|Patail200?il : hereafter 
Paper I) was dedicated to the general characteristics of LEs, 
with particular attention to the effects of self absorption on 
LE luminosity, spectral and colour appearance. In this and 
subsequent papers we will instead present and discuss appli¬ 
cations to known and test cases. In principle, LEs can affect 
all types of SNe. Due to the young stellar population from 
which they are supposed to be originated, core-collapse SNe 
are of course favored with respect to thermonuclear explo¬ 
sions. However, both due to their high homogeneity (which 
makes the analysis much simpler) and their wide use in cos¬ 
mology, here we will concentrate on Type la only. In partic¬ 
ular, the present work deals with the effects of LEs on their 
luminosity and spectra, in order both to help the observers 
in disentangling intrinsic features from LE contamination 
and to understand whether LEs can affect the conclusions 
reached on their nature. 

The paper is organized as follows. After giving a short 
introduction to the basic concepts of LE modeling in Sec¬ 
tion 1^ in Section |3 we describe the spectral library used 
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in our calculations. The results are then presented and dis¬ 
cussed in Section 2] for two different dust geometries, while 
Section is devoted to the implications on the extinction 
law. In Section 1^ we discuss our findings which are briefly 
summarized in Section [7| 


2 BASIC CONCEPTS 


The LE phenomenon has been presented and discussed in 
a number of publications and we refer to them the inter¬ 
ested reader (see for example the recent review bv ISueermarJ 
1200,?! and the references therein). The detailed description of 
the formation of a LE spectrum has been outlined in Paper 
I and, therefore, here we will only shortly recap the basic 
concepts. 

In single scattering approximation, the LE spectrum at 
a given time t can be expressed as 


Sle{X, t) 



di' 


( 1 ) 


where S°{X, t) is the intrinsic source (SN) spectrum and 
/(A, t) is a function that describes the physical and geomet¬ 
rical properties of the dust. If we define rd{X) as the optical 
depth to the SN along the line of sight and S{X, t) is the SN 
spectrum clean of LE contamination, then 


S°{X,t) = S{X,t) 


If we indicate with Cext{X) the dust extinction cross 
section, w(A) the dust albedo and g{X) the forward scattering 
degree, in all those cases when there is one single cloud with 
dust on the line of sight (like in the two examples of Fig.0, 
the formal solution of equation 0 can be written as follows 
(see Paper I): 

Sle{X, t) = G[g{X),t] LoiX) Ce.t(A) 5(A, t) (2) 


where G[( 7 (A),t] is a time- and wavelength-dependent 
function that includes dust and dust geometry properties 
and 


S{X,t) 



t') dt' 


( 3 ) 


is the time integrated SN spectrum, where tc is the time 
taken by the generic iso-delay surface to completely cross the 
dust cloud. The crossing time depends of course on the dust 
geometry and, if tc > t, the lower integration boundary of 
Equationl^has to be replaced with 0. 

When the SN is only mildly reddened {jd <d), the re¬ 
sulting LE spectrum is therefore proportional to the product 
between S{X,t) and the extinction cross section. When the 
optical depth is larger {rd > 0.3), single scattering descrip¬ 
tion must be replaced by the more realistic single scattering 
plus attenuation approximation (hereafter SSA; see Paper I), 
which holds for r 1. In this more refined schema, which 
takes into account the LE self-absorption, equation ^ be¬ 
comes 


SLEiX,t) = G[g{X),t] uj{X) Ge.t(X) 5(A,t)(4) 

where Teff{X,t) can be regarded as a weighted LE op¬ 
tical depth at a given time t. In general, one has that 
Td 7 ^ Teff, as in the case sketched in Fig. 0 (left panel). 


to the observer to the observer 



Case A 


Case B 


Figure 1. Schematic representation of dust and LE geometry for 
example cases A (large dust cloud, tc >2 yr; left panel) and B 
(circumstellar dust shell, =0.015 yr; right panel) for a generic 
observation time t {t lj2.0 yrs for Case A and t JSO.05 yrs for Case 
B). The parabolas represent the iso-delay surface corresponding 
to some preceding epochs. 


where the SN suffers a larger extinction than the LE. In 
the cases where the SN and the LE are affected roughly by 
the same extinction, the two effects tend to compensate and 
the LE spectrum is again proportional to the product be¬ 
tween the observed time integrated SN spectrum 5 and the 
extinction cross section G^xt- These formal expressions are 
useful to understand the principles of the problem and will 
be used later on to expla in the apparent peculiarity reported 
bv ISchmidt et all (119941 in the extinction function (see Sec¬ 
tion 0 . Nevertheless, we note that the exact SSA numerical 
solution is easy to implement in a code and it offers the 
advantage of properly taking into account the attenuation 
effect and the wavelength dependency of all involved quan¬ 
tities. While we refer the reader to Paper I for more details, 
here we will spend some words to discuss one important ge¬ 
ometrical effect, which is related to forward scattering and 
was only marginally mentioned in our previous paper. 

In fact, all LE calculations include the so called scat¬ 
tering phase function ‘I>(6(), which describes the probabil¬ 
ity distribi^on_of photon scattering angle 9 (see for exam¬ 
ple |c^^ahe^^8^)AJsualI^j_&is^mmtion is parameterized, 
followinB lHeiw^ fc GreensteinI (Il944l . through the forward 
scattering degree g{X). To illustrate the behaviour of $(6^), 
in Figure 0 we have plotted the results on e obtain s using 
the Milky Way dust mixtures modeled by iDraind (l2003ll 
for Rv=3.1 (solid lines) and Rv=5.5 (dashed lines) and for 
three values of the scattering angle. As one can see, the 
wavelength dependency is much stronger for lower values of 
9, with the scattering being more efficient in the blue than 
in the red. Therefore, when the dust is located far in front of 
the SN, which implies a small average scattering angle, the 
LE spectrum is expected to be bluer than in the case where 
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Figure 2. The Henyey-Greenstein scattering phase function <I>(6) 
computed for three different scattering angles and for _Ry=3.1 
and__Rv_=5^. The forward scattering degree function g is from 
|Drain3|200^ 

the same dust, with the same geometry, is placed close to 
the SN. As a consequence, the same dust cloud can produce 
the same extinction to the SN (i.e. the same rd) but different 
LE colours, according to its position along the line of sight. 
The simulations show that the [B — V) colour of the LE can 
change by ~0.15 mag from a configuration where the SN is 
embedded in the dust cloud to one where the same cloud is 
placed at more than 200 lyr in front of the SN. 


3 SPECTRAL LIBRARY 

The LE spectral modeling relies on the availability of in¬ 
put SN spectra. For this purpose one needs to build a spec¬ 
tral library with a good phase coverage and time sampling, 
especially around maximum light. This can be understood 
inspecting the cumulative functions for the template light 
curves we discussed in Paper 1. These were constructed us¬ 
ing the data of SNe 1992A and 1994 p, two standard, well 
studied and low reddening Type la iKirshner et al.l Il990t 
IPatat et al For B and V passbands, where the LE is 

expected to be brighter, it turns out that about 90% of the 
radiation is emitted in 0.2 yr, while the radiation emitted 
in the first 10 days accounts for about 5% of total. There¬ 
fore, the spectral library should include spectra from about a 
week before maximum to 2-3 months past maximum. Spec¬ 
tra at later epochs are needed only if one is going to extend 
the calculations to correspondingly late phases for dust ge¬ 
ometries with t > tc- 

We have extended the library us ed in Paper I, whic h 
included only data from SN 1998aq iBranch et ahll^OOSll . 
adding early and late spectra from SN 1992A (ESQ K ey- 
Programme on SNe) and SN 1994D llPatat et alJll99^ . so 
that the data set ranges from —11 up to -1-406 days. Irrespec- 


Figure 3. Late time spectra of SN 1992A (+328, +406) and 
SN 1994D (+293). The spectrum at +372 is an average of two 
spectra of SN 1994D and 1992A at very similar phases. The solid 
curves in the first three spectra trace the multi-Gaussian fittings. 
For the +406 spectrum, the solid line is the best fit to the +372 
spectrum scaled to match the observations. 


live of the original flux calibration, the input spectra have 
been re-calibrated to make their synthetic V magnitude con¬ 
sistent with the V template light curve. As for the late time 
spectra, we have used the available data of SNe 1992A and 
1994D, which are presented in Figure El and reach 406 days 
past B maximum light. As one can see, the spectra of the 
two objects remain very similar, also at these late stages. 
In order to remove the noise in the data, we have used a 
multi-Gaussian fit to the spectra (solid smooth lines). 

Due to the lack of data at later phases, in our calcu¬ 
lations we have assumed that after phase +406 there is no 
evolution and, from that phase on, the SN spectrum is sim¬ 
ply scaled according to the exponential decay of the V light 
curve. As a matter of fact, to our knowledge the only Type 
la SNe with optical spectroscopy at phases later than 450 
days are SNe 1991T and 1998bu, both dominated by evident 
LEs. Therefore, very little is known about the optical spec¬ 
tral behaviour at these phases. Even though there are evi¬ 
dences of deviations from the ra dioactive decay in the n ear- 
IR, as in the cases of SN 1998 bu ij Spvrom ilio et all2004^ and 
SN 2002cx llSollerman et alJl2004ll . for the optical domain 
things appear to be different. In fac t, photometric observa - 
tions of the normal la SNe 1972E iKirshner fc Okj jlIj75|L 
1992A iCaPDellaro et aP 1200111 and 1996X lls^^^^dl 
l2001^ at phases as late as 500 days indicate that the light 
curve follows the radioactive decay at least until these 
epochs. 

Another fact that should be mentioned is that even 
objects which were very similar at early phases, do 
show differences at late phases. This is clearly illus¬ 
trated in Figure [I] where we have compared the spec¬ 
tra of four Type la at about one year past maximum 
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Figure 4. Comparison between four normal Type la at late 
phases; SN1992A (ESO-Key Programme Ar chive), SN1994D 
(ESO -Key Program me Archive), SN19 98bu ISovromilio et all 
l2Q04f) and SN2001el llMattila et al.ll20Q5l) . Epochs are computed 
from the B maximum light. The thin curve in the third panel 
from bottom is a synthetic LE spectrum (see text). SNe 1998bu 
and 2001el spectra were corrected for redd ening using a ca noni- 
cal Rv=3. 1 extinction law with Ay=0.96 Jjha et al.lll999l) and 
Ay=0.57 iKrisciunas et aTll2Q03ll . respectively. 


Figure 5. Light curve calculated for a perpendicular sheet with 
a thickness AR=50 lyr placed in front of the SN with D=200 
lyr and njf=5 cm“^ {r(i{V)=0.12, solid line) and D=800 lyr and 
cm~^ (Td(\T)=0.99, long-dashed line). The short-dashed 
line traces the exponential decay of the tem plate light curve, the 
circles are the observed data for SN 1991T iSchmidt et al ]|1994 
ISnarks et al.l[T999li and filled dots mark the phases 1.2, 1.4 and 
1.6 yr. The insert shows the deviation of the measured magnitude 
from the radioactive decay tail for the first geometry. 


light: 1992A fAmisf_B)=1.47z b0.05. ^h^ips_etalJ [19991) . 
1994D rAmi^rB')=1.32 ±0.n5. IPatat et alJ Il99(j) . 1998bu 
fAmi»^f_B~l = l.nizb0.n5. [SuntzeffTl^iT 7id 2001el 

(Ami^('_B') = 1.13±0.04. iKrisciuna^^L l20o!ll . Despite the 


fact that all these SNe have shown ver 
at maximum lipht iKirshner et al.lll993l 

V similar spectra 

Patat et alill99fil 

iHernandez et alJ l200Cl. IWane et al.l l200f 

), at late phases 


some differences are visible. The most pronounced one is 
the variable intensity of the bump at a bout 7300 A, usuall y 
attributed to a blend of [Fell] lines llBowers et alJlMTill . 
which is almost invisible in SN 1992A and well developed 
in SN 2001el. Other discrepancies are observed also in the 
blue, where the [Fell] and [Felll] bumps appear to show 
different ratios in the various objects. In this respect, it is 
worthwhile noting that the spectrum of SN 1998bu is af¬ 
fected by a LE fsee lCaPDellaro et a Ll l200li : ISpvromilio' et alJ 
l2004^ , which dominates below 4500A (see Figure third 
panel, thin curve). 

The existence of these differences must be kept in mind, 
since the intrinsic properties might be mistaken with the 
onset of a LE, even when this is not the case. 


4 EFFECTS ON SN SPECTRA AND LIGHT 
CURVES 

Using the library we have presented here and the SSA ap¬ 
proximation (or the Monte Carlo approach when required) 
discussed in Paper I, we have computed the LE spectra and 
analyzed the effects these have on objects like SNe 1992A 


and 1994D under two very different conditions, i.e. when 
the dust cloud is very large (tc ^ t) and in the opposite case 
when its size is very small (tc ^ t). 

4.1 Case A: large clouds {tc ^ t) 

To illustrate the LE effects in the large cloud case (see Fig.0 
left panel), we run a series of calculations using the SSA 
approximation and a simple dust geometry, placing a dust 
sheet perpendicular to the line of sight at D=200 lyr in front 
of the SN, with a thickness AR=50 lyr and a particle density 
nrr=5 cm“®. This configuration, which implies Td{V)—0.12 
(or Av=0.13) produces a light curve which is qualitatively 
similar to what wa s observed in 1991T IlSchmidt et a l]ll994l: 
ISnarks et al.lll99'9h . as shown in Figure We remark that 
we did not make any special attempt to ht the SN 1991T 
data since our purpose here is to show the general effects in a 
realistic situation. In this respect, it is also worth mentioning 
that the light curve of this slow declining object and its spec¬ 
tra at maximum where rather different from those we have 
used for our calculations (see for example the broad light 
curve peak shown by 1991T in comparison to our template 
in Figure]^. A more detailed analysis of the SNe 1991T and 
1998bu cases is going to be presented in a separate paper. 

As shown in Figure 0 the light curve starts to deviate 
from the ordinary exponential decay (dashed line) at about 
one year after the explosion, to settle down on a very slow 
decline, totally driven by the LE, at about two years past 
the explosion. Between these two extremes, the observed 
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Figure 6. Evolution of the global spectrum (SN+LE, thick 
curve) at three different epochs for the model described in the 
text. The dashed line traces the LE spectrum while the dotted 
line is the SN spectrum reddened by Ay=0.13. For presentation, 
SN and global spectra have been normalized to their maxima. 

spectrum is expected to show a gradual transition from that 
typical of a la at late phases to that of a pure LE. 

The amount of the LE contamination can be judged in 
Figure]^ where we have plotted the results for t=1.2, 1.4 
and 1.6 yr from the explosion. Already at the first epoch, 
when the total V luminosity (SN+LE) is only 0.2 mags 
brighter than that expected for the radioactive decay, the 
resulting spectrum is sensibly different from that of a la at 
these phases, especially in the blue. Due to the peak present 
in the LE spectrum at about 4600A, the most prominent 
feature of the intrinsic nebular SN spectrum at 4700A, nor¬ 
mally attributed to [Felll], appears to be broadened and the 
region between 4000 and 4500 A much bluer than in an un¬ 
affected la. Th is is actually quite sim ilar to what was seen 
in SN 1998bu llSovromilio et af]l2004^ in the spectrum ob¬ 
tained 381 days past B maximum (see also Figure 3] here), 
where the bump at about 4000 A is unusually bright. 

The consequences of a LE in the red are milder, and 
tend to emerge only at later phases. The dominant [Felll] 
bump gets gradually fainter compared to the other emerging 
features and the colour becomes bluer and bluer, reaching 
(73 — y) ~ —0.2 when the spectrum is totally dominated by 
the LE, around 2 yr after the explosion. During the tran¬ 
sition phase a distinguishing feature is also the appearance 
of Call H&K lines, reflecting their strong presence in maxi¬ 
mum light spectra. 

The outcome can be different if the dust optical depth 
is increased, mainly due to the LE self absorption (see Paper 
I). To illustrate this fact, we have run a series of calculations 
placing the dust sheet at D=800 lyr and increasing its opti¬ 
cal depth to Tdiy) to about 1 (n_tf=40 cm“®). The distance 
was properly tuned to produce a LE with the same observed 
magnitude as the one discussed earlier (see Figure m long¬ 



time Since Explosion (yrs) 


Figure 7. Light curve calculated with the MC code for the r~^ 
wind with Ro=0-01 lyr. The curves trace the global (SN+LE, 
thick) and extincted SN (thin) light curves for Td(V)=0.5 (solid) 
and Td{V)=3.0 (dashed). The insert shows the deviation of the 
global light curve (SN+LE) from the SN light curve for the two 
optical depths. As a guide for the eye, the vertical dashed-dotted 
line marks the position of V maximum. 


dashed line). In this case the LE spectrum is much redder, 
having {B — V) ~+0.15. As a matter of fact, both the SN 
and the LE undergo roughly the same reddening, which in 
this case is E{B — V) ~0.35. 

In general, for higher extinctions the SSA approxima¬ 
tion becomes insufficient and a more complicated multiple 
scattering approach has to be performed (see Paper I). 

Before moving to the next section we like to make a 
consideration on the late spectrum shown by SN 2001el. Be¬ 
sides the intense bump at about 7500A, the region between 
4200 and 5400A appears to be rather blue. In particular, 
the bump at 4400A is much more intense than in any of the 
other objects presented in Figure |31 The synthetic colour 
of the observed spectrum is {B — y)=0.25, and it reduces 
to (B — y)=0.07 after applying the extinction reported by 
iKrisciunas et all (l2003l) for i7v=3.1. This value is unusually 
blue when compared to the ones derived from the spectra of 
SN 1994D (0.38) and SN 1992A (0.34). 

At a first glance this might be interpreted as the signa¬ 
ture of an underlying LE, but there are two quite important 
aspects that need to be considered. The first is that the 
synthetic magnitude deduced from the observed spectrum 
is y=21.24, which is in perfect agreement with the stan¬ 
dard luminosity decline. Secondly, the spectrum drops quite 
fast below 4300A, contrarily to what happens for instance to 
1998bu (see Figure |3J. The simulations show that changing 
the optical depth of the region which is responsible for the 
hypothetical LE does not allow to produce the drop which 
is seen in the blue. The conclusion is that the late time spec¬ 
tral appearance of SN 2001el is not caused by a LE of the 
type considered in this section. 
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Figure 8. Same as Figure|Hlfor Td(V)=3.0 (A^=3.26). 

In general, when the dust region is sufficiently extended 
in the direction perpendicular to the line of sight, the dust 
system integrates the emitted SN spectra over the whole 
time range from the explosion to the epoch of the obser¬ 
vation (tc is t), so that the resulting LE spectrum is very 
different from that of the SN at any given epoch. There¬ 
fore, it is very easy to detect its emergence, that can be 
difficultly mistaken with intrinsic features. An increase in 
the dust density would simply anticipate the time when the 
LE emerges, change the colour of the LE itself and cause a 
larger extinction on the SN itself, but the LE characteristic 
spectral features will still be very well pronounced. 

4.2 Case B: small clouds (tc <C t) 

More subtle is the case when the dust system character¬ 
istic dimension is significantly smaller than ct (see Eig. 0 
right panel). In fact, under these conditions, the dust sys¬ 
tem integrates the SN input spectra over times which are 
much shorter than the time elapsed from the explosion (i.e. 
tc <C t) and the resulting LE spectrum is much closer to the 
SN spectrum, making it very difficult to detect it on the basis 
of the spectral appearance alone. Of course, in order to have 
this effect, the dimension of the dust cloud must be indeed 
small. For example, for a spherical dust cloud with radius R 
and centered on the SN it is easy to show that tc = 2i?/c. 
The simulations show that in order to have a LE spectrum 
which is very similar to that of the SN at all epochs, one 
needs to have tc ^0.2 lyr, which implies R ^0.1 lyr 
cm)^. Since these dimensions are much smaller than those 
of typical interstellar dust clouds, the most natural scenario 
to be considered is the one of circumstellar material. As 

^ If one is considering only phases as late as one year, due to 
the slow spectral evolution the cloud size can be larger and still 
produce a LE spectrum very similar to that of the SN. 
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Figure 9. Evolution of the global spectrum (SN+LE, thick 
curve) at four different epochs for the wind model with iJo=0-01 
lyr and T,j(y)=0.5. The dashed line traces the LE spectrum while 
the dotted line is the SN spectrum reddened by A^=0.54. For 
comparison, the latter spectrum has also been scaled to match 
the integrated flux of the global spectrum (SN-I-LE, thin solid 
line). 

done by IChevalie'rl Jl98fJ) , IPatad (1200, and IWand i 200 ,'tl , 
we will consider here the simplest case, i.e. that of a spher¬ 
ically symmetric stellar wind in which the density is given 
by n(r) = no{Ro/r)^, where Ro is the inner boundary of the 
wind, within which the density is supposed to be null. Also, 
we assume that the wind has an outer boundary, defined 
by Ri > Rq. Under these circumstances, the optical depth 
along the line of sight is simply Td = Cext no Ro (l — Ro/Ri), 
from which it is clear that if i?i S> Ro the value of is prac¬ 
tically independent from Ri. As we have shown in Paper I 
(cf. Fig. 6 ) this geometry, compared to the others we had 
explored, is able to produce quite bright LEs with a fast 
declining light curve. 

To illustrate the effects of such a dust geometry, we 
have run a series of calculations using the MC code we have 
described in Paper F. In all simulations we have u sed the 
canon ical Milky Way dust mixture with Rv=3.1 dPraiuel 
I2n0,'d) and we have computed the LE B and V light curves 
and optical spectra for i?o=0.1, 0.01 and 0.005 lyr, keep¬ 
ing i?i = 1.0 lyr and for different values of Td{V) (0.1, 0.5, 
1.0, 2.0 and 3.0). Two examples of the resulting curves are 
presented in Fig. 0 for J?o=0.01 lyr and Td{V)=0.5, 3.0. 
The main photometric parameters obtained from the global 
light curves (SN-I-LE) are summarized in Tabled where be¬ 
sides the the optical depth and the implied extinction and 
colour excess (rd(U), Ay and E^{B — V)) we show the usual 
Amio{B), the slope at phases later than 6 months ( 7 ( 1 ?)), 

^ Proper multiple scattering treatment is mandatory in all cases 
where the dust is very close to the SN, especially when its optical 
depth is larger than 1 (see Paper I) 
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Table 1. LE effects on a dust distribution for three different values of the inner wind boundary Rq (see text). All values are in mag, 
with the exception of 'y{B) (mag/lOO*^) and Atmax (days). The number under the labels of columns 4, 5 and 6 indicate the reference 
values for the input light curve. The meaning of each column is explained in the footnotes. 


rd(V) 


E°(B - V) 

Ami5(B) 

[1.45] 

7(S) 

[1.40] 

Am3oo(B) 

[6.98] 

Atmax AVmax 

A{B - V)max 

E(B — V)max 

Ay 

A' 

AMy 







Ro = 

0.1 lyr 






0.1 

0.11 

0.04 

1.36 

1.29 

6.70 

0.0 

-0.02 

-0.01 

0.03 

0.09 

0.08 

+0.01 

0.5 

0.54 

0.18 

1.16 

1.10 

5.91 

1.2 

-0.10 

-0.05 

0.13 

0.44 

0.40 

+0.04 

1.0 

1.09 

0.35 

0.95 

1.01 

5.23 

2.6 

-0.24 

-0.07 

0.29 

0.84 

0.88 

-0.04 

2.0 

2.17 

0.70 

0.55 

0.93 

4.29 

3.4 

-0.48 

-0.16 

0.54 

1.69 

1.68 

+0.01 

3.0 

3.26 

1.05 

0.30 

0.89 

3.65 

4.8 

-0.75 

-0.25 

0.80 

2.52 

2.48 

+0.04 







i?o=0.01 lyr 






0.1 

0.11 

0.04 

1.38 

1.41 

6.94 

0.0 

-0.05 

-0.02 

0.02 

0.05 

0.05 

0.00 

0.5 

0.54 

0.18 

1.30 

1.43 

6.79 

1.2 

-0.27 

-0.09 

0.09 

0.28 

0.28 

0.00 

1.0 

1.09 

0.35 

1.24 

1.43 

6.63 

2.6 

-0.56 

-0.14 

0.21 

0.53 

0.64 

-0.11 

2.0 

2.17 

0.70 

0.97 

1.37 

6.31 

3.2 

-1.07 

-0.28 

0.42 

1.10 

1.29 

-0.19 

3.0 

3.26 

1.05 

0.86 

1.33 

6.04 

4.4 

-1.56 

-0.41 

0.64 

1.70 

1.99 

-0.29 







Ho=0.005 lyr 






0.1 

0.11 

0.04 

1.41 

1.41 

6.96 

0.0 

-0.06 

-0.02 

0.02 

0.00 

0.05 

0.00 

0.5 

0.54 

0.18 

1.39 

1.45 

6.88 

1.1 

-0.30 

-0.09 

0.08 

0.24 

0.25 

-0.01 

1.0 

1.09 

0.35 

1.29 

1.48 

6.78 

1.5 

-0.60 

-0.17 

0.18 

0.49 

0.55 

-0.06 

2.0 

2.17 

0.70 

1.19 

1.49 

6.60 

2.8 

-1.17 

-0.30 

0.40 

1.00 

1.23 

-0.23 

3.0 

3.26 

1.05 

1.03 

1.47 

6.46 

3.1 

-1.69 

-0.44 

0.61 

1.57 

1.90 

-0.33 


A^ = 1.086Td(V'); E°{B - V')=A^/3.1; Ami5(_B) = Bl5 - Bmax', 'y{B) = {dB/dt)t^lSO', Am 30 o(-B) = S 300 - B^ax', ^tmax=tmax - Rnax'^ 
AV^ax=V^ax - A(B - VUax=(B - V)max - {B - + E°(B - V); E(B - V)=E°(B - V) + A{B - VUax', 

Av=A°y + AVruax-, A'y=3AE{B - V); AMv=Av - A'y. 


the magnitude jump between maximum and 300 days past 
maximum {Amzoo(B)), the time shift in the B maximum 
{Atmax), the differences in V and {B — V) at maximum 
with respect to the extincted SN {AVmax, A{B — V)max), 
the colour excess with respect to the unreddened SN at 
maximum {E{B — V)max), the apparent extinction obtained 
comparing with the intrinsic magnitude {Av), the extinction 
derived from E{B — V)max assuming Rv=2>A (Ay) and the 
deviation of the absolute magnitude at maximum with re¬ 
spect to the input value (AMy). For comparison, the cor¬ 
responding values for the template light curves assumed by 
the models are as follows: Ami^{B)=lAb mag, '^(B)=1AQ 
mag/lOO"^ and Am3oo(B)=6.98 mag. 

There are several conclusions that can be drawn already 
inspectin g Table [J Fir st of all, as already qualitatively no¬ 
ticed by IWa.ngi feOPfil . the resulting light curves become 
broader and the maximum is attained later then in the pure 
SN case (cf. also Fig. |7|l. This has the clear effect of de¬ 
creasing Ami 5 (i 3 ) rather sensibly and the change increases 
with Td{V) and i?o. Moreover, the global colour gets bluer 
than that of the purely extincted SN and this effect becomes 
stronger for larger Td and smaller Rq. In principle, unless the 
values of this parameters become too deviant from th e nor¬ 
mal ones (0.9 ^ Ami 5 (B) <1.9. I^illios et alJll999ll . they 
could be attributed to an intrinsic feature of the SN. Never¬ 
theless, there are other side effects that tend to betray the 
presence of an underlying LE. In fact, the late time decay 
rate 7 (B) increases with Td and Ro, while the magnitude 
jump Am 3 oo(B) grows with Rq and decreases for increas¬ 
ing values of Td- In the case of 7?o=0.005 lyr (^5x10^® cm). 


the global light curve around maximum light changes but 
its overall shape remains similar to the template light curve 
(e.g. giving plausible values for Ami^(B), Amzoo(B) and 
so on). This has a very interesting consequence, at least 
for Ro SiO.Ol lyr. In fact, besides measuring a slower post¬ 
maximum decline rate, one would also measure a colour 
which is significantly bluer and a global apparent luminos¬ 
ity which is higher than those of the purely reddened SN. 
As the simulations show, these two facts combined together 
would lead an hypothetical observer deducing the reddening 
correction from the global observed colour (and using a stan¬ 
dard value for Ry) to derive an absolute luminosity which 
is brighter than the real one by an amount that can reach 
about 0.3 mag (see last column of Tab. 0. In conclusion, one 
would attribute the smaller value of Am\^(B) to a larger in- 
trins ic luminosity, as in t he case of the classical Phillips rela¬ 
tion iPhillins et al.ll99flll . Of course this does not mean that 
the observed relation is produced by dust in the vicinity of 
the SN, since Amis correlates in fact wi t h other features like 
spectral line ratios ||NugenL.et alJ[l9i^ |Benettj,.eL_^| 20 ()j: 

iHachinger, Mazzali fc BenettiTT^OOnBong^^^^l .1 l200.»jl . 

But certainly, the presence of hidden LEs might introduce an 
additional spread in the measured photometric parameters. 

The effects on spectra are presented in Fig. and 
where we show the time evolution in the Bo= 0.01 lyr case 
for r£i(P)0=0.5 and 3.0. In the low optical depth case the 
effect is rather mild (the global spectrum is slightly bluer 
than that of the reddened SN spectrum at all phases), while 
the LE influence becomes much more relevant in the high 
optical depth case. As a consequence of the increased LE 
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contribution to the total flux, the time integration becomes 
visible and this causes the spectral features to show clear 
deviations from the intrinsic ones. The most remarkable is 
the absorption profile of the Si II 6355 A line around max¬ 
imum light (t=0.05 lyr), which is much broader. The same 
behaviour is shown by all lines and at all epochs. A late 
phases, the line peaks also appear blue-shifted with respect 
to the intrinsic features. Moreover, as in the low density case, 
the global spectrum is bluer than the reddened SN spectrum 
at all times. 

In general, the strongest effect is expected around max¬ 
imum light and in the weeks that immediately follow it, i.e. 
when the spectral evolution is fast compared with the LE 
integration time which, in the case of i?o=0.01 lyr is of the 
order of 12 days®. In the pre-maximum phase, the early spec¬ 
tra contribute quite marginally to the global flux, due to the 
very rapid rise in luminosity while, just after maximum, the 
SN flux in the preceding days is actually higher than on the 
current epoch, so that the contribution is indeed relevant. 
This is clearly displayed in Fig. EOl where we present the 
calculated evolution of the Si II region between t=0.03 and 
t=0.10 yrs, which correspond to —7.3 and -1-18.3 days past 
maximum light. On the first epoch, the profile of the global 
spectrum is already broader than in the input SN, but the 
deviation grows in the following days, reaching a maximum 
during the two weeks past maximum light. The position of 
the absorption through minimum appears to be red-shifted 
with respect to the input spectrum, producing a lower pho- 
tospheric velocity estimate. During the second week, also 
the emission peak of the line is affected and it appears much 
broader, reflecting the higher velocities reached in the pre¬ 
ceding epochs. These deviations, together with those shown 
by the light (and colour) curves, would most probably allow 
one to recognize the presence of a LE. Nevertheless, we must 
notice that this is true only for high optical depths, while 
for rd{V) ^1, the effects on photometric parameters and 
spectral features are rather mild (see Tableland Fig. EJ, 
making the detection of the LE certainly much more diffi¬ 
cult. In those cases the net effect of the LE would be to alter 
the colours and the overall flux of t he obs e rved spectrum, 
with the consequences discussed bv IWan3 ||200^, that we 
will analyze in the next section. 


5 IMPLICATIONS ON EXTINCTION AND 
EXTINCTION LAW 

If a SN is affected by the emergence of a LE, it is clear that 
each time one tries to derive the amount of extinction or 
information about the dust properties, the results will be in 
general misleading. As we have done in the previous section, 
we distinguish here the two opposite cases. 


5.1 Case A 

When l^chmidt et alJ (Il994l discovered the LE in SN 1991T, 
they noticed that, in order to reproduce the observations, 

® Due to the r~^ dust density distribution, this is larger than the 
7.3 days one would have in the case all the dust was located in a 
very thin shell with a radius of 0.01 lyr. 


they had to correct the time integrated SN spectrum S for 
the effects of scattering. These were parameterized with a 
function of the type A““ and the best fit was obtained using 
a=2. As we have seen in Sectional the correction one has 
to apply to S includes three components: the first is related 
to the light scattered by the dust into the observer’s line 
of sight (directly proportional to Cext), the second accounts 
for the reddening suffered by the SN (proportional to e*^®**) 
while the third describes the auto-absorption within the dust 
cloud that generates the LE itself (proportional to 
Finally, there is a mild wavelength dependency introduced 
by the albedo and a more important contribution by the 
function G{X,t) (see SectionEJ. A more convenient formula¬ 
tion can be obtained passing from the extinction cross sec¬ 
tion Cext{X) to the normalized extinction law a(A). If we 
pose 


a(A) 


CextiX) 

axt(v) 


then the absorption (in magnitudes) can be expressed 
as A(A) = a(X)Av and equation ||1J can be rewritten as 

e ' (5) 

where k(V) = Ad{V) — Aeff{V) is expressed in magni¬ 
tudes. In a first approximation, one can parameterize a(A) 
with a power law of the type X~^. From equation 0 it is 
clear that a fitting to the ratio between Sle and S using a 
power law A““ would give a value for a which is not directly 
related to /3, since other terms with different wavelength de¬ 
pendencies enter its expression. In fact, using our simula¬ 
tions for the perpendicular sheet with D=200 lyr we have 
found that fo r the Milky Way dust mixture with i7v=3.1 
the best fit values of a depend mostly on 
the value of k{V), while the dependence on r^ff is very 
mild, at least in the validity range of the SSA approxima¬ 
tion (r^ff ^1). Typical values for a are 1.4, 1.7, 2.0 and 2.2 
for k( 1/)=—0.2, 0.0, 0.2 and 0.4 respectively. 

It is interesting to notice that the a=2 value found 
by ISchmidt et alJ (ll994^ can be explained with a canoni¬ 
cal extinction law and a small difference between the ex¬ 
tinction suffered by the LE an d the SN, i.e. onfy 0 .2 mag. 
The SN was certainly reddened: IPhillins et alJ I^1999^ report 
in fact a total extinction Av=0.53±0.17, 0.46 of which is 
due to the host galaxy. If we assume Ad=0A6 (or equiv¬ 
alently Td{V)=QA2) we can therefore give an estimate to 
Aeff, which turns out to be 0.26 (i.e. re//(E)=0.24). The 
corresponding model is presented in Figure mi Given the 
fact that the synthetic spectrum has been computed us¬ 
ing the spectra of objects like 1994D and 1998aq which are 
quite different from SN 1991T, especially during the max¬ 
imum light phase, the match is reasonably good. However, 
we must notice that the observed spectrum is also compati¬ 
ble with k{V)=0, i.e. with the SN and the LE suffering the 
same amount of reddening fFigure ITTl thin smooth curve). 

Another ex ample of a possibly simila r situation is given 
by SN 1998bu. ICappellaro et afT J200lll have presented a 
LE spectrum at 670 days past maximum light, with a 
{B — V) colour close to 0, i.e. very s imilar to that observed 
for SN 1991T iSchmidt et alJfl994^ . On the other hand, 
l.lha et aP (|l99^ have reported a tot al extinction Av=0.94 
±0.15, while iHernandez et alJ (l2000h found Ay=1.0 ±0.1. 
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Rest Wavelength (A) 

Figure 10. Evolution of the Si II 6355A region for ilo=0.01 lyr 
and Tjj(y)=3.0 between —7.3 to +18.3 days past maximum light. 
The dotted curves trace the reddened SN spectra. For presenta¬ 
tion they have been scaled by a suitable amount in order to match 
the global spectrum. The vertical dashed line marks the position 
of the absorption minimum of Si II in the first spectrum. In the 
MC simulations the photons where collected with a resolution of 
AA=15 A and a time bin AT=0.01 yrs. 

Using the same geometrical configuration, the best fit with 
the observed LE spectrum is obtained for Te//(U)=0.62, 
which implies Av=0-67 (see Fignre I14II . Given the fact 
that the estimated galactic extinction in the direction of 
SN 1 998bu is only Ay=0.08 llSchlegel. Finkbeiner fc Davis! 
^9^, this value is about 0.2 mag lower than the average 
value reported by IjhT et al.l 1 I 1999 II and iHernandez et all 

J200(t . 

In this respect it is interesting to note that, due to the 
geometrical effect related to forward scattering that we have 
described in Section |21 if the dust is close to the SN, then 
a lower dust optical depth would still produce a red LE 
spectrum. This is not because there is more reddening, but 
because under these circumstances the scattering efficiency 
does not grow at shorter wavelengths as fast as it does if 
the dust is located far way. This is illustrated in the insert 
of Fig.im where we have plotted the results of a simulation 
where the same perpendicular sheet is placed at distance 
D—Q. In this case the best fit is obtained for ffiV)=0.25, 
which would imply a SN extinction Av~0.27, which is sig¬ 
nificantly lower than the values reported in the literature. 

In conclusion, the information on the extinction and the 
extinction law one derives from the observed LE spectrum is 
geometry dependent and not univocally related to the dust 
properties. 

5.2 Case B 

As we have seen in Sec. 10 the presence of a LE produced 
by a dust geometry with tc is able to introduce sensible 



Rest Wavelength (A) 

Figure 11. Spectrum of SN 1991T at 729 days past B maxi¬ 
mum (ESO Key Programme). The s pectrum was correct ed for 
the Milky Way reddening 1A^T=0■07. iPhillios et ^ll999^ ■ The 
thick smooth curve is a synthetic LE spectrum obtained with the 
spectral library described in the text with Te//=0.24. For com¬ 
parison, a synthetic spectrum with r^f f=0A2 is also shown (thin 
curve). 


luminosity and colour changes with mild effects on the other 
photometric and spe ctroscopic cha racteristics, provided that 
'^d{V) ^ 1- This, as IWan j has shown, can have an 

interesting consequence on the properties one would deduce 
for the dust. In fact, if we take the example of iilo=0.01 lyr 
and Td{V)=l.Q, from Table^we see that, due to the presence 
of a LE, the global V magnitude is 0.56 brighter than the 
purely extincted SN, which is affected by an extinction of 
Ay—1.09 mag. Therefore, assuming a distance to the SN and 
using the expected absolute luminosity, one would deduce 
Ay— 0.53. Then, comparing the colour at maximum with a 
reference value (which here we assume to be the one of our 
template object), one would also deduce E{B— V)max=0.21 
from which she would erroneously conclude that Rv ~2.5. 
As far as the photometric properties are concerned, such 
an object would show normal Ami 5 {B) (1.24) and 7(5) 
(1.43) and at late phases the slight over-luminosity (0.3 mag) 
would probably go unnoticed due to the measurement errors. 
Therefore, at least on the basis of those parameters, it would 
be very difficult to recognize the presence of an underlying 
LE. 

Since the difference between the global luminosity and 
that of the pure SN is not constant in time and changes 
rather rapidly around maximum light (see Fig.|71 upper right 
insert), the values of Ry deduced in this way are expected 
to change with time in a similar way. This is illustrated 
in Fig. where we present the results of MC simulations 
for 5o=0.01 lyr and two different values of the dust optical 
depth. The total-to-selective absorption ratio changes in fact 
very rapidly from the explosion to about 0.15 yrs, which cor¬ 
respond to the maximum of the intrinsic SN {B — V) colour 
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Figure 12. Upper panel: colour evolution for the wind dust ge¬ 
ometry with Ro=0.01 lyr for t^{V)=1.0 (lower thick curve) and 
T^(y)=3.0 (upper thick curve). The dashed lines trace the corre¬ 
sponding colour curve in the case of pure extinction and i?ix=3.1, 
while the dotted curve is the SN intrinsic colour. Lower panel: ra¬ 
tio between Av = ^E{B — V)max deduced comparing 
the global colour with the intrinsic colour of the SN at maximum 
for T^(V)=1.0 (solid line) and 3.0 (dashed line). The insert shows 
the behaviour of Ay in the two cases. In all plots the vertical 
dotted line marks the epoch of maximum light. 

curve (dotted curve). After that epoch, the ratio remains 
almost constant. This behavio ur is qualitatively similar to 
the one found bv lWan3 (l200fill (see his Fig. 2) and the dis¬ 
crepancies are probably due to the different treatments of 
multiple scattering^. 

It is important to remark that the global colour curve is 
not obtained by a simple rigid reddening of the unextincted 
SN colour curve. On the contrary, the difference between 
bluest and reddest colour is smaller and the tail at phases 
later than 0.15 yrs (~55 days) is less steep. The effect be¬ 
comes more pronounced for higher optical depths. In other 
words, once the observed global colour curve is corrected for 
the colour excess deduced at maximum, this is not going 
to match a reference colour curve and the late behaviour 
would clear ly differ from the homogeneous trend found by 
iLiral 1^9^. We reckon that this col our eff e ct, sh own also 
by the semi-analytical calculation of IWanel l^2005^ (see his 
Fig. 3, lower panel) is a good diagnostic to unveil the pres¬ 
ence of a hidden LE. 

At this point we must emphasize that, in all practical 
cases, the colour excess is deduced from the data using the 
colour curve. On the contrary, the extinction is usually de¬ 
rived from this value after assuming Rv and then applied to 
the observed magnitude to o btain the absolute lu minosity of 
the object (see for example iBenetti et al.lEnn4l . Following 

^ We reach results very similar to those reported bylWan^ {200^ 
for Td(V')=0.5. 


Figure 13. Upper panels: comparison between the global spec¬ 
trum (thick curve) and the intrinsic SN spectrum (thin curve). 
Lower panels: comparison between intrinsic (i7ix=3.1, dashed 
curve) and derived (solid curve) normalized extinction laws. The 
simulations were performed using I?o=0.01 lyr and Tci{V)=1.0. 


this procedure in the case of our previous example would 
lead to a result which is actually not far from the real one. 
In fact, if one adopts E{B — y)=0.21 and assumes Rv=3.1, 
this would produce an absolute magnitude which is only 0.11 
mag brighter than the real value (see Tabled last column). 
Remarkably, this is of the sam e order of the correct ion fore¬ 
seen by the Phillips relation llPhillips et al.lll999ll for the 
given Ami 5 (i?). As we had already noticed in Sec. 14.11 this 
effect is most likely going to add a spread to the observed 
luminosity-decline rate relation. 

Another way of deriving information on the extinction 
law A(A) is to compare the observed spectrum S{X) with 
a reference unreddened SN spectrum 5° (A) obtained at the 
same phase and corrected to the same distance (see for in¬ 
stance lElias-Rosa et aljl200,‘j) . If the LE effects were negli¬ 
gible, one could in fact write A{\) = 1.086 ln[S'°(A)/S'(A)] 
and hence derive directly the extinction wavelength depen¬ 
dency. Of course, when ^(A) contains a significant contribu¬ 
tion from the LE the result is going to be misleading. This 
is clearly illustrated in Fig. US where we present the results 
obtained from our Monte Carlo simulations at two differ¬ 
ent epochs (t=0.05 and 0.15 yrs) for the case t(V)=1.0 and 
7?o=0.01 lyr. For the first epoch, which approximately cor¬ 
responds to the B maximum light, the colour of the global 
spectrum is clearly redder than that of the unreddened SN 
(E(B — V) ~0.2, see also Fig. ll2ll and the derived extinction 
law is just slightly steeper than the input one, thus giving a 
smaller Rv- With the exception of the Sill 6355A feature, 
the derived curve is reasonably smooth. Totally different is 
the case at t=0.15, which roughly corresponds to the max¬ 
imum of the {B — V) colour curve. At this phase, in fact, 
the colours of the global spectrum and the unreddened SN 
are similar (see also Fig. ll2ll and the derived extinction law 
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is dominated by the differences displayed by the two spec¬ 
tra due to the time integration effects. Of course, in a real 
case, this could be interpreted as the presence of a LE but 
also with a deviation from the standard spectral evolution. 
Irrespective of the physical reason, the apparent evolution 
shown by the derived extinction law is a clear diagnostic 
signaling that additional effects are present and that Ry 
can not be estimated with this method. Given these facts, 
a detailed analysis of extincted la should be able to clarify 
whether this mechanism is indeed at work and whether it 
can explain the low values of Rv derived by several authors. 


6 DISCUSSION 

All known cases of LEs in Type la SNe have revealed them¬ 
selves through the gradual appearance of an additive flux, 
during the nebular phase, in the spectral region bluewards 
of 4500A, where one can easily recognize the imprints of the 
time integrated SN spectrum. Therefore, they are instances 
of the Case A we have described in Sec. 14.11 The spectrum 
discussed bv ISovromilio et al.l (l2004^ is a very good exam¬ 
ple of this behaviour (see also Figure |1] here). Moreover, due 
to the close wavelength coincidence of the dominating [Fe 
III] nebular feature with the most prominent bump of a LE 
spectrum, the former line is broadened. Besides this, a dis¬ 
tinguishing spectral feature of a LE is the clear presence of 
the Ca II H&K lines which, though, at the beginning of the 
transition phase can be masked by the presence of a similar 
structure observed in pure nebular spectra at about 4000 A. 

In general it is fairly easy to recognize the Case A LE 
contamination on the basis of the late spectra alone, espe¬ 
cially if the optical depth is not too high, so that its spectrum 
is definitely blue, with (B — U) ^0. In cases where the LE 
forms in regions of high dust density, its spectrum can turn 
red without necessarily having the SN spectrum reddened. 
In this case, the LE contamination is a bit more subtle and 
tends to leave unchanged the overall colour of the observed 
spectrum. However, the characteristic LE bumps at about 
5900A and 6500A (see Figure]^ appear quite clearly. 

An interesting fact that emerges from the analysis of 
the two known cases of SNe 1991T and 1998bu when com¬ 
pared to our synthetic LE spectra is that clear differences 
are visible. As far as 1991T is concerned, its LE dominated 
spectrum taken at 729 days past maximum (cfr. Figure rTTIl 
shows interesting features, which are directly related to its 
peculi ar appearance during th e maximum phase (see for in¬ 
stance In particular, we notice that 

the absence of the distinguishing la feature (i.e. the Si II 
6150 A) during the pre maximum phase and its weakness in 
post maximum is reflected in the LE spectrum. In fact, while 
our synthetic spectra, calculated using “normal” events like 
SNe 1994D and 1998aq, do show a clear imprint of the P-Cyg 
profile displayed by this line during the photospheric phase, 
this is practically absent in the data of 1991T. Another clear 
difference is that visible in correspondence of the Call H&K 
lines. Again, the weakness of this feature in 1991T is a conse¬ 
quence of what happened in the pre-maximum pha se, when 
this line was totally absent dFilippenko et af]ll992^ . 

In the other case, the match between the synthetic and 
the observed LE spectrum is fairly good but, still, the ob¬ 
served deviations can all be attributed to small differences 


in the spectral appearance of SN 19 98bu with respect to 
SN 1994D fsee iHernandez et alJl200(1i . 

All of this has an interesting consequence, i.e. that the 
maximum phase peculiarities can be easily seen in the LE 
spectrum. Therefore, not only one could classify an historical 
SN on the basis of its LE spectrum alone (disentangling for 
example between a Type I and a Type H), but one could also 
get some more detailed information about its features during 
maximum light, so that it would be possible to distinguish, 
for instance, between a normal la, a 1991T-like or a 1991bg- 
like event. This might have some relevance when traces of an 
historical SN explosion are found by detecting its LEs, as in 
the re cent cases nea r the SN remnant of Cas A llKra.iise et alJ 
l2005^ and in LMC iR.est et al.ll2005)l . Provided that the LE 
is bright enough to allow a spectroscopic observation, this 
wo uld provide a powerf ul too l to type the SN. T he detections 
by iKrause et al.l (l2005^ and iRest et ^ (l2005^ lets us hope 
that more cases will be discovered in connection with SN 
remnants in our Galaxy, offering a unique chance to directly 
classify the explosion that generated them and to establish 
a firm connection between the two. 

In principle one could envisage to use the Case A LE 
observations to derive some information about the extinc¬ 
tion law in the host galaxy, simply comparing the SN time 
integrated spectrum and the LE spectrum. In fact, under 
the conservative assumption that the LE and the SN suf¬ 
fer the same extinction (i.e. k(U)=0 in equation 0), the 
ratio between Sle and S is reddening independent. Unfor¬ 
tunately, even in these simplified circumstances, one is left 
with the albedo and forward scattering degree wavelength 
dependence. More in general, due to the way the effects of ge¬ 
ometry, optical depth and dust properties combine together, 
one cannot derive unambiguous information on the extinc¬ 
tion law since very similar solutions can be obtained us¬ 
ing different ingredients. Some meaningful results might be 
achievable from spectroscopic observations of LEs if some of 
the dependencies can be removed by other and independent 
observations. Of great help, in this respect, are high spatial 
resolution images, which can give direct information on the 
geometrical dust distribution. 

With the increased interest on Type la SNe related to 
their usage in cosmology and the availability of larger and 
larger observing facilities, it is most likely that an increasing 
number of objects belonging to this class will be observed 
in the near future and followed up to even more advanced 
phases than it is currently done. Our forecast is that these 
investigations will probably be disturbed by the presence of 
LEs, too faint to be detected by the typical current instru¬ 
mentation. In fact, as we have shown here, a dust density 
of the order of n_ff=5 cm“® is sufficient to produce a LE 
compatible with those seen in SNe 1991T and 1998bu. Due 
to the lack of LE detections in the vast majority of the cases 
jBoffi et al in Paper I we had concluded that the 

height scale of Type la must be significantly higher than that 
of the dust in the disk of h ost galaxies. Indeed, a very sim ilar 
conclusion was reached bv IPella Valle fc Panaei^ (ll992^ on 
the basis of completely independent considerations. 

Nevertheless, when one pushes the observations at later 
phases, due to the fast decay typical of these objects during 
the nebular phase (~5 mag yr“^), the light curve rapidly 
reaches levels where a LE produced by a much lower den¬ 
sities would dominate the spectrum. For example, for the 
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Figure 14. Spectrum of S N 1998bu at 670 days past B maximum iCannellaro et fd1l2nmh . The spectrum was corrected for the Milky 
Way reddening fv4T/=0.08. Ijha et ‘ailllQQflh . The smooth curves are synthetic LE spectra with f =0.62 (thick) and Tf,ff=0.90 (thin) 
respectively. The upper right insert shows the results for the same dust sheet placed at D=0 and Tf,f f=0.25 (Ay=0.27). 


same geometry used in Figure^ a density of n_ff=0.5 cm“® 
would produce a LE approximately 2.5 mag fainter, imply¬ 
ing that the transition to a LE dominated spectrum would 
have started only 6 months later. Therefore, the study of 
Type la at phases later than 1.5-2 yr, aimed for example 
at detecting deviations from the radioactive decay, will be 
most probably hindered by the emergence of a Case A LE. 

If, on the one side, it is true that there is an increasing 
number of instances where LEs of this kind are actually ob- 
served (see also the recent case reported bv ISugermanlboO^ 


and IVan Dvk. Li fc Filiopenk^ l2005l for the core-collapse 
SN 2003gd), one may in fact wonder whether and to which 
extent early phases could be affected, modifying for instance 
the light curve shape, the overall luminosity and introducing 
peculiarities in the spectral appearance. As we have shown 
in this paper, this can actually take place when the dust 
cloud is very close to the SN and has small dimensions, cir¬ 
cumstances that we have generically in dicated as Ca se B. 
As it has been independently shown bv IWan3 (|200^, cir- 
cumstellar dust can in principle alter the observed light and 


























































































































Reflections on Reflexions 13 


colour curves producing, for example, interesting effects on 
the derived extinction curve. Our analysis (see Secs. 14. H and 
15.211 . run under the simplifying hypothesis that the dust is 
distributed according to a law with an inner cavity with 
radius Ro, shows that already with moderate optical depths 
{Td{V)=l), the presence of dust in the immediate vicinity of 
a Type la {Ro cm) would turn into significant changes 

in the photometric parameters, which would still be within 
the normally observed ranges. In particular, the post maxi¬ 
mum decline rate Amis would be smaller and the SN would 
also appear inherently brighter, due to the additional flux 
contribution added by the underlying LE, i.e. producing an 
additional scatter in the Pskovski-Phillips relation. In this 
respect, we notice that there are many instances of Type la 
showing low decline rates (Ami 5 (B) <1) without displaying 
the spectral peculiarities ty pical of slow-declinin g SN 1991T- 
like events (see for example IPhillins et aI.llT999l . Table 2). 

Besides influencing the photometric parameters, the 
time integration during the fast evolving phase around max¬ 
imum light tends to produce changes in the line widths and 
peaks, altering the measured expansion velocities. As we 
have said, according to our simulations the deviations from 
a LE-free standard behaviour would be rather subtle and 
could be easily misidentified as the intrinsic features of a 
peculiar object. Whilst most of the photometric parameters 
tend to be mildly altered (at least for relatively low values 
of the dust optical depth) the colour evolution shows quite a 
clear deviation from the homogeneous behaviour show n by 
Type la between 30 and 90 days past maximum light IILiral 
|l99^. If indeed all LE-free objects follow that path, then 
the colour curye could be a good diagnostic tool to unyeil 
the presence of a LE. _ _ 


As pointed out by IWand ||200^, the presence of 
Case B LEs could explain the low yalues for Rv which 
haye been deduced from SN observations (see for example 


PhilliDS et alJl999llKriscmnas et alJ200r 

■iKnoD et alJ200.Sl 

Altavilla et alJl2004l andlKrisciunas et al 

120051 for the most 


recent examples). As we have shown here, this is in prin¬ 
ciple true, even though some side effects should be clearly 
detectable, at least for the simple dust distributions we have 
analyzed. 

According t o the l o wer p anel of Fig. 1121 or to the equiv¬ 
alent Fig. 3 of IWangl fePOdf . one would expect that the 
deduced Rv remains roughly constant (~1) in this phase 
range. Nevertheless, in those plots E{B — V) is kept at 
the constant value attained at maximum light, while the 
difference in colour between the global spectrum and the 
SN spectrum changes with time, both in our Monte Carlo 
calculations and in the semi-analytical model presented by 
IWangl { 200 ^. In particular, when the global colour becomes 
equal to or even bluer than that of the unreddened SN, the 
extinction correction is expected to produce weird results. 
Therefore, if one attributes the deviations from the standard 
colours to a different extinction law, she also has to accept 
the puzzling conclusion that this has to change with time. 
We reckon that both the colour evolution and the time vari¬ 
ability of the derived extinction law would be sufficient to 
infer that the observed deviations are not due to a peculiar 
dust composition but rather to a LE. 

Another interesting fact worth to be mentioned is the 
discrepancy that one would And between the extinction de¬ 
rived from the equivalent width of the Nal D lines and 


that deduced from photometric or spectroscopic compar¬ 
isons with templates. In fact, while the intensity of the Nal D 
absorption is related only to the dust present along the line 
of sight (and would hence give a direct estimate of r^), all 
other methods would be affected by the LE contribution. As 
we have seen in Sec. E21 this tends to significantly decrease 
the amount of derived extinction. In the case of normal dust- 
to-gas ratios and for a simple r~^ spherically symmetric 
dust geometry, this would imply that the values deduced 
from the Nal doublet are systematically higher than those 
inferred from other methods (compare Ay with Av or Ay 
in Table Ql. 

In all this discussion we have assumed (and required) 
that there can be dust very close to the SN. In this respect 
there is a very important consideration one has to make. 
In fact, if on the one hand it is true that LE effects be¬ 
come more subtle when the dust is closer to the SN (and 
can therefore be mistaken with intrinsic properties), on the 
other there actually might be a lower limit for the SN-dust 
distance. In fact, for typical expansion velocity of 10"* km 
s“^, the ejecta would reach i?o=10^® cm in less than four 
months. This is a problem that needs to be addressed with 
detailed modeling, in order to see whether the interaction 
would be able to destroy and/or sweep away the dust and 
to evaluate possible effects on the observed reddening, its 
time evolution and so on. One should also mention that the 
radiation field produced by the SN explosion would probably 
evaporate the dust up to a certain radius. If in the case of 
a core-collapse event the UV flash is expected to produce a 
cavity in the du st with a radius of cm (see for example 

IChevali^IlQSfifl . the lack of UV photons in a Type la might 
still make the dust survive the explosion. Another implica¬ 
tion that needs to be addressed is related to the fact that 
substantial amounts of dust should be associated with cor¬ 
respondingly high amounts of gas. If the gas is close enough 
to the SN, the n some interaction is expected, like in the case 
of SN 2002ic JHamuv et al.ll200^ . 

Finally, what remains to be studied is the effect of asym¬ 
metric dust distributions (disks, tori, lobes and so on) that, 
among other things, would certainly produce a net contin¬ 
uum polarization up to several per cent (see Paper I). Be¬ 
ing close to the SN the polarized flux could be in fact high 
enough to produce detectable effects and serve as a useful di¬ 
agnostic tool to identify the presence of an underlying Case 
B LE. Finally, to make the scenario fully self-consistent, one 
should also investigate the effects of dust heating on IR re¬ 
emission, which might turn out to be another observable 
side-effect produced by this kind of LEs. 

Detailed analysis of reddened objects will certainly clar¬ 
ify whether Case B LEs are indeed taking place and this, in 
turn, will probably give us a few more hints about the na¬ 
ture of Type la progenitors. We are currently workin g on the 
highly reddened SN 2003cg llElias-R.osa et a^]l2005^ , which 
shows clear features that could be interpreted as the signa¬ 
ture of a Case B LE. The results will be reported in another 
paper of this series. 


7 CONCLUSIONS 

In this paper we have analyzed and discussed the effects of 
LEs on light curves, colours and spectra of Type la SNe for 
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two different dust geometries: extended clouds (Case A) and 
small circumstellar clouds (Case B). The main results can 
be summarized as follows: 

(i) The effect of Case A LEs is limited to the late phases. 
The exact time when the LE contribution dominates the 
global luminosity depends on dust distance and density. 

(ii) There is a transition phase, during which the LEs con¬ 
tribution becomes more and more important and the global 
spectrum turns from a normal nebular spectrum into a pure 
LE spectrum. 

(iii) During the transition phase, the overall colour tends 
to become bluer, spectral features broader and emission line 
peaks displaced with respect to a pure SN nebular spectrum. 

(iv) During the LE-dominated phase the spectrum retains 
the peculiarities shown by the SN at maximum, so that the 
LE detection in an historical SN allows to classify it with a 
good confidence level. 

(v) The spectroscopic observation of Case A LEs in the 
vicinity of SN remnants can help in establishisng a firm and 
direct relation between the SN type and the remnant. 

(vi) The information on the extinction properties one can 
derive from LEs (both Case A and B) observations may lead 
to incorrect conclusions on the dust nature. 

(vii) Case B LEs affect in principle all phases, without 
necessarily producing notable effects on the spectra. This 
is particularly true if the dust is confined within ~10^® cm 
from the SN explosion and Td <1. 

(viii) In the presence of circumstellar dust, the light 
curves tend to become broader with a smaller Amis. The 
rising time increases, the late time decay becomes slower 
and the luminosity difference between maximum light and 
nebular phase tends to decrease. 

(ix) The SN absolute luminosity deduced from the appar¬ 
ent magnitude, colour excess and a standard extinction law 
is brighter than the real value. 

(x) The difference in luminosity is of the same order of the 
correction foreseen by the Pskowski-Phillips relation. This, 
coupled to the fact that Amis is decreased, is a possible 
source of noise in th e luminosity- decline rate relation. 

(xi) As shown by IWand (1200,di , Case B LEs modify the 
luminosity and colours in such a way that the apparent 
selective-to-total extinction ratio Rv is significantly de¬ 
creased. 

(xii) In this scenario, Rv should show a rapid evolution 
during the first months after explosion and this, together 
with a non standard colour evolution, should be a simple 
diagnostic tool to unveil the underlying LE. 
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